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Abstract 


An appropriate flow field in the bipolar plates of a fuel cell can effectively enhance the reactant transport rates and liquid water removal efficiency, 
improving cell performance. This paper proposes a novel serpentine-baffle flow field (SBFF) design to improve the cell performance compared 
to that for a conventional serpentine flow field (SFF). A three-dimensional model is used to analyze the reactant and product transport and the 
electrochemical reactions in the cell. The results show that at high operating voltages, the conventional design and the baffled design have the 
same performance, because the electrochemical rate is low and only a small amount of oxygen is consumed, so the oxygen transport rates for 
both designs are sufficient to maintain the reaction rates. However, at low operating voltages, the baffled design shows better performance than 
the conventional design. Analyses of the local transport phenomena in the cell indicate that the baffled design induces larger pressure differences 
between adjacent flow channels over the entire electrode surface than does the conventional design, enhancing under-rib convection through the 
electrode porous layer. The under-rib convection increases the mass transport rates of the reactants and products to and from the catalyst layer and 
reduces the amount of liquid water trapped in the porous electrode. The baffled design increases the limiting current density and improves the cell 


performance relative to conventional design. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


Fuel cells are electrochemical reactors now used for a wide 
variety of applications. Due to their high efficiency (nearly twice 
that of the present generation of internal combustion engines), 
portability and near-zero emissions, fuel cells are attractive as a 
power source for automobiles. Of the various types of fuel cells, 
the proton exchange membrane fuel cell (PEMFC) operates at 
near-room temperatures and is considered to be a good choice 
for automotive applications. However, the commercialization of 
the PEMFC is still hindered by several technological problems, 
among which is severe water flooding of the cathode and the 
resulting mass transport losses [1—3]. Over recent decades, in an 
effort to improve PEMFC performance, many analyses, models 
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and numerical simulations have been developed [1-20] for vari- 
ous transport phenomena and electrochemical kinetics to gain a 
better understanding and to develop strategies for optimal design 
and operation scenarios. 

The flow field design in the bipolar plates is one of the most 
important issues in a PEMFC. An appropriate flow field design 
in the bipolar plates can improve the reactant transport and the 
efficiency of the thermal and water management. To this end, 
different flow field configurations, including parallel, serpentine, 
interdigitated, and many other combined versions, have been 
developed. Many efforts have been devoted to the optimal flow 
field design to improve cell performance [14-33]. 

Nguyen [24] first proposed the interdigitated flow field, in 
which baffles were added at the ends of some channels. The 
baffle design forced the reactants to flow through the gas dif- 
fusion layer (GDL) and the shear force of this reactants flow 
helped blow out the liquid water trapped in the inner lay- 
ers of the electrodes. As a result, the mass transport rates of 
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Nomenclature 


z>] 


SETTAT 


gE STEND DOI 


ac] 


chemical activity of water vapor 

cathode cross-sectional inlet flow area (m2) 
reference exchange current density (A m7?) 
reaction area (m?) 

source term of variable é 

mass fraction 

quadratic drag factor 

hydraulic diameter of the flow channel (m) 
mass diffusivity (m? s7!) 

Faraday constant (96,487 C mol— ty 

current density (A m~?) 

current density (A m7?) 

coefficient of water vapor condensation rate (s7 1) 
coefficient of water vapor evaporation rate (s7!) 
permeability (m7) 

distance along the channel measured form the 
inlet (m) 

molecular weight (kg mol~!) 

pressure (atm) 

mass flow rate (kg s7! m7?) 

universal gas constant (8.314 J mol`! K7!) 
volume ratio occupied by liquid water 

time (s) 

temperature (K) 

X-direction velocity (m s7!) 

Y-direction velocity (m s7!) 

fuel velocity at the cathode inlet (m s7!) 
Z-direction velocity (m s7!) 

cathode pressure drop loss (W m7?) 

width of rib (m) 

X-direction coordinate (m) 

Y-direction coordinate (m) 

Z-direction coordinates (m) 

species valence 


Greek letters 


Oa electrical transfer coefficient in forward reaction 
Oe electrical transfer coefficient in backward reaction 
E porosity 

n overpotential (V) 

À water content in membrane 

u viscosity (Pa s) 

v kinematic viscosity (m? s7!) 

ay exchange coefficient 

p density (kg m~?) 

Om electric conductivity of membrane 

T tortuosity of the pores in the porous medium 

0) dependent variables 

® phase potential function (V) 

superscripts 

ref reference value 


subscripts 

a anode 

aver average 

c cathode 

eff effective 

in inlet 

k kth species of the mixture 


H+ hydrogen ion 
H2 hydrogen 
H20 water 

O2 oxygen 

sat saturation 
total total 


the reactants from the flow channel to the inner catalyst layer 
were improved and the water flooding problem at the cath- 
ode was significantly reduced. Soong et al. [31] proposed a 
novel flow channel configuration in which the baffles were 
inserted transversely in the channel of a conventional parallel 
flow field to form a partially blocked fuel channel. They found 
that reducing the gap size and/or increasing the number of baf- 
fles enhanced the reactant transport and cell performance. Liu 
et al. [32] investigated the effect of baffle-blocked channels on 
the reactant transport and cell performance of PEMFC with a 
conventional parallel flow field. Their results indicated that the 
baffles forced more reactants from the flow channel into the GDL 
to enhance the chemical reactions which augmented the cell 
performance. 

These investigations have shown that enhancing the con- 
vection of reactants through the GDL by adding baffles is 
an effective way to reduce water flooding at the cathode and 
increase the reactant mass transport, thus improving both the 
cell performance and the operating stability. Based on this 
understanding, the effects of three new serpentine-baffle flow 
field designs for the PEMFC were investigated in this study. 
The flows were obtained by inserting the baffles in different 
locations in the channels for a conventional serpentine flow 
field. A three-dimensional numerical model was used to com- 
pare the cell performance for baffled and conventional designs. 
The effects of liquid water formation on the reactant trans- 
port were taken into account in the model. The oxygen mass 
flow rates and the liquid water distributions at the interface 
between the cathode GDL and the catalyst layer (CL), and 
local current densities for baffled and conventional designs are 
analyzed to show the advantages of the baffled design. The com- 
pressor powers for baffled and conventional designs are also 
evaluated. 


2. Flow field design 


Water flooding normally occurs at the cathode electrode of 
the PEMFC because the electrochemical reaction on the cathode 
produces water vapor. If the partial pressure of the water vapor 
is higher than the saturation pressure, water vapor condenses to 
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Fig. 1. Cathode flow fields designs. (a) SFF; (b) SBFF-1; (c) SBFF-2; and (d) SBFF-3. 


form liquid water. When a large amount of liquid water accumu- 
lates in the porous layer pores, the oxygen transport resistance 
increases and the oxygen mass flow rate decreases. Therefore, 
the cathode flow field design is a key factor for enhancing reac- 
tant and product transport and for removing the liquid water. The 
present paper analyzes various cathode flow field designs. The 
results are compared to a triple serpentine flow field to identify 
the advantages of the novel baffled design. 

Fig. 1 compares the cathode flow fields for the conventional 
and baffled designs. Fig. 1(a) shows a conventional triple ser- 
pentine flow field, which includes three serpentine loops, AA’, 
BB’, and CC’. The reactants enter the cell from the inlets of 
the three loops and then flow along the channels. Some of the 
reactants diffuse into the GDL and CL and are consumed by 
the electrochemical reactions, while the remainder flows out of 
the cell from the outlets of the three loops. Fig. 1(b)—-(d) shows 
three triple serpentine-baffle flow field designs, named SBFF- 
1, SBFF-2, and SBFF-3. The SBFF-1 includes three serpentine 
loops, DD’, EE’, and FF’, with the same inlet configuration as 
the conventional design, but with two baffles added at the outlets 


of DD’ and FF’. Therefore, the reactants only flow out of the cell 
from the outlet of EE’. The SBFF-2 design has three serpentine 
loops, GG’, HH’, and II’, with three baffles added at the inlet of 
HH’, and at the outlets of GG’ and II’. The SBFF-3 design also 
has three serpentine loops, JJ’, KK’, and LL’, with two baffles 
added at the inlet and outlet of KK’, and another two baffles 
placed at the center of JJ’ and LL’ as shown in Fig. 1(d). The 
anodes for all four cell designs are parallel flow field designs 
with 12 flow channels and 11 ribs. 

This paper considered miniature fuel cells with dimensions 
of 23 mm x 23 mm x 2.845 mm. All four designs had the same 
23 mm x 23 mm reaction area dimensions, and the same GDL, 
CL, and membrane thicknesses. The detailed physical dimen- 
sions of the fuel cells are summarized in Table 1. The operating 
conditions were all the same for a fair comparison for all the 
cells. The fuel cell temperature was assumed to be 323 K, the 
reactants on the anode side included hydrogen and water vapor 
with a relative humidity of 100%, the reactants on the cathode 
side contained oxygen, nitrogen, and water vapor with a rel- 
ative humidity of 100%, the inlet flow rate on the anode side 
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Table 1 

Fuel cell dimensions 

Quantity Value 
Reaction area (mm?) 23 x 23 
Channel width (mm) 1 

Rib width (mm) Í 
Channel height (mm) 1 

Rib height (mm) 1 
Anode GDL thickness (mm) 0.4 
Anode CL thickness (mm) 0.005 
Membrane thickness (mm) 0.035 
Cathode GDL thickness (mm) 0.005 
Cathode CL thickness (mm) 0.4 


was 260 cm? min™!, the inlet flow rate on the cathode side was 


700 cm? min™!. 


3. Model 


A three-dimensional model of the full fuel cell was used to 
analyze the electrochemical reactions and transport phenomena 
of the reactants and products. The cell was divided into the anode 
flow channels, membrane electrode assembly (MEA, including 
the anode GDL, anode CL, proton exchange membrane, cathode 
CL, and cathode GDL), and cathode flow channels. The gov- 
erning equations included the mass, momentum, species, and 
electrical potential conservation equations. The model assumed 
that the system is three-dimensional and steady, the inlet reac- 
tants are ideal gases, the system is isothermal, the flow is laminar, 
the fluid is incompressible, the thermophysical properties are 
constant, and the porous GDL, CL, and proton exchange mem- 
brane (PEM) layers are isotropic. 

The transport equations for the three-dimensional PEMFCs 
are given as 

Continuum equation: 


ðu ðv dw 
a a A i 1 
TT E “ 


Momentum equation: 


Eeff OP u u Pu 
Sys ae ( poj ) +S (2) 
IX IX Z 


p a əx? ay?" az? 
i v Es aon ov 
u— — +w— 
cae Gey a 
Eeff OP vu v dv 
F ay + ear (5 + ay2 + a2 + Sy ( ) 


Eeff OP w w dw 
=~ + veer (Ta + - ) +5. 4) 


Species equation: 


ac, ac aC, 
cor (u ate Faw =) 


x dy az 
aCe aCe 3Cr 
=D Se + Si 5 
cer ( az t my + ae ) + Sc + SL (5) 


In the momentum equation, € is the porosity and S,, Sy, and 
Sw the corrected terms of the reactant flow in the gas diffusion 
layer and the catalyst layer and of the proton transfer in the 
PEM, which is listed in Table 2. In the species equation, Dj. cfr 
is the effective diffusion coefficient and Se represents the source 
terms due to the chemical reaction in the catalyst layer and the 
proton exchange membrane. The Bruggeman correction [34] 
is employed to describe the influence of the porosity on the 
diffusion coefficient: 


Dyett = De” (6) 


If the partial pressure of water vapor is greater than the satura- 
tion pressure, the liquid water forms in the PEMFC. To consider 
the effect of the liquid water formation, itis assumed that the pore 
in the porous material is blocked by liquid water, which results in 
the modification of the diffusion coefficient and the porosity in 
the species equation. The source term Sy, in the species equation 
due to the liquid water is determined by [35] 


€eftCH,0 : 
y= Mmoke— pr PO — Psat), if Puyo > Psat (7) 


keEetts( Psat — P50); if PHO < Prat 


where the saturation s is the ratio of the liquid water volume to 
pore volume in the porous material, M the molecular weight of 
water, ke the condensation rate constant of water, ke the evapo- 
ration rate constant of water, and £eff is the modified effective 
porosity of the porous medium by considering the liquid water 
effect which is given by 


Eeff = e(1 — s) (8) 
The saturation pressure of water can be expressed as 

= = -572 —IT3 
Poat = 1072 17940.02953T—9.1837x 107% T?-+1.4454x 10-77 (9) 


To calculate the local current density, the phase potential 
equations should be solved 


a a® if ə a® P a a® 5 (10) 
or or or = S; 
ax | T ax ay \ T ay az | az j 


where S;=0 in the membrane, $;=—ja in the anode catalyst 
layer, S; =je in the cathode catalyst layer, ® the phase potential, 
and om is the ionic conductivity of the membrane. In this study, 
Butler—Volmer equation [36] is used to calculate the transfer 
current density generated by the electrochemical reaction: 


Cy 1 

~ ya ef | “Ho (OaF/RT)n _ __ 

Ja = Ajoa (S) f JET. a) 
2 


Cc 1 
- 4 tef | Or (@aF/RT)y oO 
Je = Ajo (=) le sare (12) 
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) 


where A get is the reference exchange current density, œa and œe 


the transfer coefficient in anode and cathode catalyst layers, 
the overpotential, R the gas constant, and T is the temperature 
of the fuel cell. The ionic conductivity of the membrane can be 


Zo , o , o 
7 


Ke calculated by [37] 
1 1 
ae, ae ae T) = o!" exp |1268 ( — — — 13 
2 E42 E Gi al 303 T aa 
z ERE ë 
a Š and the reference ionic conductivity is 
oe on (®) 
ol >o È ô $ SÈ otf = 0.005139\ — 0.00326 (14) 
4 — J 0.043 + 17.814 — 39.85a° + 36.007 O<a<l 15) 
~ ) 1441.4(a-1) l<a<3 
where a is the activity of water vapor at the cathode side. The 
relationship between the phase potential and the current density 
+ is given by 
. o@ 
y ly = mg (16) 
Š a 
Pp 
ee iy = —Om— (17) 
3 ti oy 
9 -e ] Ñ 
3 P a Se op 
i E iz = —Om— (18) 
dy 


Boundary conditions at the anode flow channels and the cath- 
ode flow channels are as follows: the inlet flow rates are constant, 
the inlet gas compositions are constant, and the flows are fully 
developed at the outlets of the anode and cathode flow channels. 
At the solid walls, no slip and zero fluxes are hold. At the inter- 
faces between the gas channels, the diffuser layers, the catalyst 
layers, and the PEM, equalities of the velocity, mass fraction, 


Š i momentum flux, and mass flux are applied. More details were 
a given elsewhere [30-33]. All parameters used in the model are 
i X 1 È listed in Table 3. 
we aie gile N Since these equations for this complex convection—diffusion 
Blo Pale problem cannot be solved analytically, it was solved using the 


finite volume method on a collocated cell-centered grid. The 


Table 3 
Parameters of the fuel cell 


Parameter Value 
Aig 9 x 108 
Aig 1.5 x 10? 
Qa 0.5 
a = Oe 1.5 
£ Ss Echannel 1 
3. ` Tchannel 1 
a S Kchannel oo 
= R N EGDL 0.5 
lalo -21> TGDL 1.5 
5, kept 1.76 x 107! 
2 ECL 0.4 
z TCL 1.5 
2 7 kcL 1.76 x 1071! 
3 z 5 EMem 0.28 
a8 z E E TMem Dagan model 
s 2 £ Q H 2 kMem 1.8 x 10718 
EA Oo © ©) = 
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general convection—diffusion equation can be expressed in con- 
servative form as 


ə = 
3 (OP) + V(pud — EgVo) = So (19) 


where ¢ is the time, ġ the general dependent variable, &4 the 
exchange coefficient, Sy the source term, ú the velocity vector, 
and p is the density. When Eq. (20) is integrated over a control 
volume, the resulting finite volume equation is 


app = degre + aww + angy + ass + andy + aot +b 
(20) 


where @p is the value of variable ¢ at node P of the control vol- 
ume, dz, . . ., ỌL the values of the variable at nodes in neighboring 
control volumes, ap, . . ., a coefficients in the discretized equa- 
tions and b is the source terms for the discretized equation for 
ġ. The semi-implicit method for pressure-linked equation (SIM- 
PLE) algorithm, developed by Patankar [38], was employed to 
solve the governing equations. The coupled set of equations was 
solved iteratively, with the solution considered to be converged 
when the relative error in each field between two consecutive 
iterations was less than 10~°. In this work, the CPU time (Intel 
Core 2 Duo E6300) for each case was about 20h. 

The model used non-uniformly distributed elements with 
116, 93, and 33 elements in the X-, Y- and Z-direction (the flow 
channel inlets are along the Y-direction and the height is in the Z- 
direction). The grid independence was examined in preliminary 
test runs. Four non-uniformly distributed grid configurations 
were evaluated for the PEMFC with the parallel flow channel 
design. The numbers of elements in the X-, Y- and Z-direction 
were: (I) 70 x 70 x 25, (ID 93 x 93 x 33, (IID) 116 x 93 x 33, 
and (IV) 116 x 116 x 41. The influence of the number of ele- 
ments on the polarization curves (J—V¢e| curves) for the fuel cell 
is shown in Fig. 2. The polarization curve for grid (I) differs 
from those for grids (II-IV). At an operating voltage of 0.3 V, 
the deviation is about 3.9%. However, the calculated polarization 
curves for grids (II), (IID and (IV) do not show any notable differ- 
ences. At an operating voltage of 0.3 V, the difference between 
the current densities for (II) and (III) is about 0.32% and the 
difference for (III) and (IV) is about 0.25%. Therefore, grid (IID) 
was chosen for the simulations as a trade off between accuracy 
and execution time. 
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--0-- 93x93x33 
ssssAses 116x93x33 
- ~~ -116x116x41 


V ceil (V) 


` 0 3000 6000 9000 
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Fig. 2. Influence of the number of elements on the polarization curves. 


4. Under-rib convection effect 


When the reactants are pumped to flow through the flow chan- 
nel from the inlet to the outlet, the local pressure in the channel 
can be approximated as [20,27-29,39] 


128uQ 
4 
ndee 


P(L) = Pin — L (21) 
where Pin is the inlet pressure of the reactants, P(L) the local 
pressure of the reactants, u the reactants viscosity, Q the mass 
flow rate of the reactants, der¢ the effective hydraulic diameter of 
the flow channel, and L is distance along the channel measured 
from the inlet. Eq. (1) indicates that the pressure decreases with 
downstream due to viscous losses. For the conventional and baf- 
fled designs, each serpentine loop proceeds through three 180° 
turns to form 12 channels and 11 ribs. For each rib, the two 
channels on its left and right sides may belong to the same loop, 
or may belong to two different loops; therefore, the pressure 
differences at each location between the two channels on the 
left and right sides differ for the various cases. For example, for 
the conventional design shown in Fig. 1(a), the two channels 
on the left and right sides of rib 1 belong to the loops AA’ and 
BB’. The inlet pressures of the reactants of these two channels 
are almost equal; thus, the pressure difference between points a 
and b is nearly zero. However, the two channels on the left and 
right sides of rib 2 belong to the same loop, CC’; therefore, the 
pressure at point c is higher than that at point d. The pressure dif- 
ference across the rib, i.e., between points c and d, is expressed 
as 


(La — Le) (22) 


Therefore, the pressure difference between points c and d is 
proportional to their distance, Ly—L,. When the two points are 
near the channel inlet, the pressure difference is maximized, but 
is minimized near the channel bend. When the pressure gradient 
across the rib is much larger than the pressure gradient along 
the channel direction, considerable cross-leakage flow of reac- 
tants from point c to point d occurs between adjacent channels 
(referred to as under-rib convection [25]). According to Darcy’s 
law, the average superficial velocity in the porous layer across 
the rib can be expressed as 


= — z (ha — Le) (23) 
T 


where kp is the permeability of the reactants through the GDL 
and W, is the rib width. Therefore, when the inlet flow rate of 
the reactants, the channel and rib dimensions, and permeability 
of the GDL remain constant, Vayer is proportional to the distance 
between points c and d, La—Le. Therefore, the under-rib con- 
vection is stronger near the inlet of loops CC’ and weaker at 
the channel bend. For the conventional design, under-rib con- 
vection also occurs under ribs 3 and 4 as well as under rib 2. 
The under-rib convection induces strong convection in the elec- 
trode, bringing the reactants to the CL and removing the liquid 
water from the reaction sites and porous electrodes, in a way 
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similar to the flow in the interdigitated flow channels, which 
results in the better overall cell performance with the serpentine 
flow channels. Therefore, this under-rib convection significantly 
influences the PEM fuel cell performance with the serpentine 
flow field design. 

The baffles change the pressure distribution of the reactants 
in the channels, which will increase under-rib convection. Fur- 
thermore, near the baffles, as in the interdigitated flow field 
design, the reactant transport is converted to forced convec- 
tion. Therefore, the baffled design is expected to enhance the 
reactant transport rates. For the SBFF-1 design, with the con- 
ventional design, the under-rib convection occurs under the ribs 
5, 6, and 7. Moreover, because the baffles are added at the outlets 
of loops DD’ and FF’, the reactants still have higher pressures 
near the outlets of DD’ and FF’. For example, the pressures at 
points g and i are greater than at point h. Thus under-rib con- 
vection also occurs under ribs 8 and 9. The SBFF-2 design also 
has a baffle at the inlet of loop HH’; therefore, the pressures 
at points j and 1 are greater than at point k, and the pressure 
at point m is greater than at point n. The changes in the pres- 
sure distributions in adjacent channels results in more under-rib 
convection at more ribs, including ribs 10, 11, 12, 13, 14, 15, 
16, and 17. Using a similar analysis for the SBFF-3, under-rib 
convection occurs at ribs 18, 19, 20, 21, 22, 23, 24, 25, and 
26. 


5. Results and discussions 


To verify baffled design does improve cell performance, the 
I-V¢ey polarization curves and J-Wee power density curves for 
the conventional and baffled designs are given in Fig. 3. For 
operating voltages greater than 0.7 V, the conventional and baf- 
fled designs have the same performance, indicating that the cell 
performance is independent of the flow channel design. For oper- 
ating voltages lower than 0.7 V, the cell performances of the 
conventional and baffled designs differ, the differences increases 
with the deceasing operating voltage. At the same operating volt- 
age, the total current densities and power densities for the three 
baffled designs are all higher than for conventional design. The 
cell performance best for the SBFF-3 design, then the SBFF-2 
design, the SBFF-1 design and finally the conventional design. 
For example, at operating voltage of 0.3 V, the total current den- 
sities for SBFF-1, SBFF-2, and SBFF-3 are 23,983, 24,882, and 
25,499 A m2, which means that the total current densities were 
increased by about 10.36%, 14.50%, and 17.34% compared with 
the 21,731 A m~? for the conventional design. Thus, the baf- 
fled designs show superior performance over the conventional 
design. 

Because the flow field designs significantly affect cell perfor- 
mance at low operating voltages, the local transport phenomena 
for the operating voltage of 0.3 V were analyzed further to under- 
stand the effect of the baffles on the cell performance. Since the 
power output is the consequence of the electrochemical reac- 
tion, the consumption of oxygen along the cathode GDL-CL 
interface is an index of the cell performance. Higher oxygen 
mass flow rates along the cathode GDL—CL interface means that 
more oxygen enters the CL to participate in the electrochemi- 
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Fig. 3. Conventional and baffled cell performance characteristics. 


cal reaction per unit time, resulting in higher current densities 
and better cell performance. The distributions of the oxygen 
mass flow rates along the cathode GDL-CL interface for the 
conventional and baffled designs are presented in Fig. 4. For 
the conventional design (Fig. 4(a)), the maximum oxygen mass 
flow rate occurs at the flow channel inlet and decreases along the 
channels forwards the outlets. Since the oxygen is gradually con- 
sumed and liquid water accumulates due to the electrochemical 
reaction along the flow channels. However, at the channel bends, 
the oxygen mass flow rate has a local peak, because the corners 
force oxygen flow into the GDL by forced convection. Fig. 4(a) 
also indicates that under ribs 2, 3, and 4, where the under-rib 
convection occurs, the oxygen mass flow rates are significantly 
higher, while under other ribs for the conventional design, the 
oxygen mass flow rates are very low. The higher pressure dif- 
ferences across ribs induce the stronger under-rib convection, 
so the local oxygen mass flow rates at these locations are larger. 
For example, for rib 4, the oxygen mass flow rates at point e 
are higher than at point f. Comparison of Fig. 4(a) and (b-d) 
shows that the oxygen mass flow rates in the baffled designs are 
remarkably different from that in the conventional design. The 
different flow rates result from two main differences between 
the baffled and conventional designs. First, near the baffles the 
oxygen transport is changed from diffusion to forced convection 
due to the baffle blockage effect; therefore, more oxygen enters 
the GDL, inducing the increased local oxygen mass flow rate. 
Secondly, the baffles change the reactant pressure distribution, 
which causes the additional under-rib convection and increases 
the local oxygen mass flow rates in the GDL under the ribs. For 
example, for SBFF-1, besides under ribs 5, 6, 7 (same locations 
as ribs 2, 3, and 4 for the SFF), the oxygen mass flow rates also 
significantly increase under ribs 8, and 9 due to the convection 
effects under the ribs when the flow field with baffle design is 
adopted. Similarly, for SBFF-2, additional the convection occurs 
under ribs 10, 11, 13, 16, and 17 (ribs 12, 14, and 15 have the 
same locations as ribs 2, 3, and 4 for the SFF), for SBFF-3, addi- 
tional the convection occurs under ribs 18, 19, 21, 23, 25, and 
26 (ribs 20, 22, 24 have the same locations as ribs 2, 3, and 4 for 
the SFF), which causes that the oxygen mass flow rates signif- 
icantly increase under these ribs. Therefore, the baffled design 
enhances the oxygen utilization efficiency compared with the 
conventional design. 
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Fig. 4. Oxygen mass flow rate (kg m7? s) distributions on the cathode GDL-CL interface. (a) SFF; (b) SBFF-1; (c) SBFF-2; (d) SBFF-3; and (e) y= 11.5 mm. 


The oxygen mass flow rate also significantly affects the 
liquid water distribution. High oxygen mass flow rates mean 
that more oxygen enters the CL per unit time to participate 
in the electrochemical reaction and more liquid water is pro- 
duced. On the other hand, high oxygen mass flow rates also 
induce large shear forces which more efficiently remove the liq- 
uid water. Fig. 5 shows the liquid water distributions on the 
cathode GDL-CL interface for the conventional and baffled 
designs. Comparison of Figs. 5 and 4 shows that the liquid water 
distributions are opposite to the oxygen mass flow rate distri- 
butions, indicating that increasing the oxygen mass flow rate 
tends to promote liquid water removal. Liquid management is 


very important for improving PEMFC performance. If the liquid 
water produced by the electrochemical reaction is not removed, 
it will accumulate in the pores of the porous layer which signif- 
icantly increases the oxygen transport resistance and decreases 
the electrochemical reaction rate. The cell performance sharply 
decreases when flooding occurs. Thus, appropriate flow field 
designs can enhance the reactant transport rate and the liquid 
water removal rate, as a result, the mass limitations are reduced 
and cell performance improves. Fig. 5 shows that the baffled 
designs have higher oxygen mass transport rates and liquid water 
removal efficiencies than the conventional design; therefore, the 
liquid water content in the porous layers is significantly reduced. 
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Fig. 5. Liquid water distributions on the cathode GDL-CL interface. (a) SFF; (b) SBFF-1; (c) SBFF-2; (d) SBFF-3; and (e) y= 11.5 mm. 


Fig. 6 presents the local current densities on the middle cross- 
section in the PEM for the conventional and baffled designs. The 
oxygen consumption rates determine the power output; thus, a 
high oxygen mass flow rate corresponds to a high local current 
density. Comparison of Figs. 6 and 4 shows that the distribu- 
tions of the local current densities are similar to those of the 
oxygen mass flow rates, confirming the accuracy of the model. 
Figs. 4—6 indicate that with the baffles, the baffled designs sig- 
nificantly improve the oxygen utilizations, by producing higher 
oxygen mass flow rates and lower liquid water distributions in 
the porous layers; therefore, the cell performance significantly 
improves. Because the baffles are added at different locations in 


the channels, the strength and extent of the additional under-rib 
convection between adjacent channels and the forced convection 
near the baffles differ for the three baffled designs, resulting in 
different cell performance characteristics for SBFF-1, SBFF-2, 
and SBFF-3. On the whole, the performance of SBFF-3 is better 
than that of SBFF-2, which is better than that of SBFF-1. 

The flow field design influences not only the cell performance 
but also the pressure drop in the fuel cell. Larger pressure drops 
in the fuel cell mean that more power is needed to pump the 
reactants. Thus, the pressure drop is a significant issue to be 
considered in choosing the flow field designs in addition to the 
I-Vecei Curve. Figs. 7 and 8 present, respectively, the pressure 
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Fig. 6. Current density distributions (A m~*) on the middle cross-section in the PEM. (a) SFF; (b) SBFF-1; (c) SBFF-2; (d) SBFF-3; and (e) y= 11.5 mm. 


distributions on the middle cross-section in the cathode flow 
channels and the total pressure drops on the cathode side for 
the conventional and baffled designs. The total pressure drops 
are, respectively, 608, 1204, 817, and 692 Pa for SFF, SBFF-1, 
SBFF-2, and SBFF-3. It indicates that the pressure drops for 
baffled designs are higher than that for the conventional design. 
The cathode pressure drop is related to the power density by 
[21]: 


A PA cha V 
Wp car (24) 
Atotal 


where Wp, represents the cathode pressure drop loss, AP 
the total cathode pressure drop across the fuel cell, Acha the 
cathode cross-sectional inlet flow area, V the fuel velocity 
at the cathode inlet, and Ajota; is the reaction area. The cal- 
culated cathode pressure drop losses and output power for 
the conventional and baffled designs are listed in Table 4. 
Although the pressure drops for the baffled designs are larger 
than for the conventional design, the pressure drop losses are 
far less than the cell output power for the miniature fuel 
cells analyzed here. Thus, the pressure drop losses can be 
ignored. 
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Fig. 8. Total pressure drops on the cathode side for the conventional and baffled 
designs. 


Table 4 
Estimation of pressure drop losses at the operating voltage of 0.3 V 


Flow field type Ap (Pa) Ween (Wm~?) Wy (Wm?) Wret (Wm?) 


SFF 608 6518.96 4.47 6515.49 
SBFF-1 1204 7194.98 8.85 7186.13 
SBFF-2 817 7464.65 9.01 7455.64 
SBFF-3 692 7649.92 7.63 7642.29 


6. Conclusions 


An appropriate design of the flow channels in the bipolar 
plates can significantly impact the thermal and water manage- 


ment of the PEMFC. A novel baffled design is presented here 
to provide the significantly better performance than the conven- 
tional design. A three-dimensional numerical model is used to 
predict the local transport phenomena and cell performance for 
the baffled and conventional designs. The conclusions drawn 
from the analyses are 


1. The analysis shows that the effect of the baffles on cell per- 
formance is dependent on the operating conditions. At high 
operating voltages, the flow field design has little effect on 
cell performance because the electrochemical reaction rate is 
low and only a limited amount of oxygen is consumed with 
only a small amount of liquid water produced. Therefore, the 
oxygen transport rates for both the baffled and conventional 
designs are sufficient to maintain the reaction rates. 

2. At low operating voltages, the cell performance is strongly 
dependent on the flow field design. Compared with the con- 
ventional design, the baffled design significantly enhances 
the mass transport capability and the efficiency of liquid water 
removal, which results in better performance than with the 
conventional design. The baffles in the channels on the cath- 
ode side force the reactants near the baffles into the GDL and 
CL to enhance the electrochemical reaction rate at the CL. 
The baffles also increase the reactant pressure in the chan- 
nels; thus, the under-rib convection occurs under more ribs, 
which increases the reactant mass flow rate under those ribs. 

3. For the miniature PEMFC analyzed here, the compressor 
power needed to overcome the cell pressure drop is far less 
than the cell output power and can be neglected; thus, the baf- 
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fled design has superior performance over the conventional 
design. 
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